Effective Onsite Interaction for Dynamical Mean-Field Theory 
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A scheme to incorporate non-local polarizations into the dynamical mean-field theory (DMFT) 
and a tailor-made way to determine the effective interaction for the DMFT are systematically in- 
vestigated. Applying it to the two-dimensional Hubbard model, we find that non-local polarizations 
induce a non-trivial filling-dependent anti-screening effect for the effective interaction. The present 
scheme combined with density functional theory offers an ab initio way to derive effective onsite in- 
teractions for the impurity problem in DMFT. We apply it to SrVOa and find that the anti-screening 
competes with the screening caused by the off-site interaction. 

PACS numbers: 71.10.-w, 71.27.+a 



-Introduction. Understanding physical properties of 
strongly correlated electron systems is one of the most 
challenging subjects in condensed matter physics 0, ■ 
For this purpose, it is essential to capture fermionic 
many-body effects necessitating a proper and accurate 
treatment of a large number of interacting fermions. The 
large number of electronic degrees of freedom in real 
materials are intractable, even with rapidly developing 
computational power. Hence, various ingenious ways of 
reducing the degrees of freedom have been developed. 
Aside from the reduction to mean-field effective one- 
particle Hamiltonians, as in density functional theory 
(DFT) , including dynamical fluctuations for the reduced 
and tractable degrees of freedom is a route that has been 
explored extensively over the last decades. 

Approaches have been proposed [H-Q to partially trace 
out the degrees of freedom far from the Fermi level, leav- 
ing an effective low-energy model for a small number of 
bands near the Fermi level. The resulting Hubbard-type 
lattice fermion models are much simpler than the origi- 
nal problem containing a huge number of bands. This re- 
duction (downfolding) has been successfully incorporated 
in the constrained random phase approximation (cRPA) 
[3| by the use of maximally localized Wannier orbitals 
(MLWO) [f| as a basis set. It should be noted that, by 
tracing out certain electronic degrees of freedom, the ef- 
fective interactions in the lattice fermion models (e.g. the 
Hubbard U, as exemplified by jj cKPA derived with the 
cRPA) are much reduced compared to the original bare 
Coulomb interactions [6-14 [ because of the screening by 
polarizations of the eliminated degrees of freedom. 

Although several efficient ways to solve the lattice 
fermion models have been proposed [lj], it is still too 
difficult to treat realistic situations so that a further re- 



duction is highly desired. The widely used dynamical 
mean- field theory (DMFT) [H, ^ indeed offers a prac- 
tical way of describing local correlation effects along this 
line [H, where the lattice fermion models are mapped 
onto quantum impurity models. 

Although U cHPA is widely used as input for DMFT 
calculations, the conventional cRPA treatment totally 
excludes non-local screening processes within the target 
band. These are also not contained in the DMFT, which 
only accounts for the local screening processes. Hence, in 
the present work we argue that a better starting point is 
the inclusion of non-local screening processes of the target 
band within the RPA yielding an effective onsite inter- 
action [/ DMFT . Albeit tailor-made interaction parameters 
for the impurity problem were employed in Rcf. Il7l a 
systematic investigation has been missing so far. 

In this Letter, we examine a scheme for the systematic 
determination of the effective onsite interaction [/ DMFT 
for DMFT calculations. This scheme is applied to both 
the two-dimensional (2D) single-band Hubbard model 
and to SrV03 by using an ab initio description. The 
application to the Hubbard model unexpectedly reveals 
the inequality (J BMFT > U and a non-trivial filling depen- 
dence of {J DMFT with a peak around the van Hove singu- 
larity. A filling-dependent f/ DMFT is also observed in the 
ab initio results for SrVC>3. These are ascribed to an 
anti-screening effect induced by non-local polarizations, 
namely, a test-charge electron induces an off-site hole or 
electron and they again induce an onsite electron. This 
nonlocal effect increases £/ DMFT . The present elucidation 
contributes not only to the specific determination of the 
DMFT-interaction parameters, but also to gain insight 
into the nature of the reduced and simplified fermionic 
models in general. 
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-Equations to derive {7 DMFT . Here, we derive the basic 
equations to evaluate [/ DMFr from first principles calcula- 
tions [17f . In the RPA, the screened Coulomb interaction 
W can be written as (1 — v Xo)~ lv with the independent- 
particle polarization xo and the bare Coulomb interac- 
tion v. The polarization xo is divided into xb an d Xo; 
where xb is a polarization formed in the target subspace 
and Xo is the res t- Note that this decomposition is not 
necessarily restricted to bands (cRPA); it is also applica- 
ble to the real space using localized basis sets. For exam- 
ple, the "dimensional downfolding" has been formulated 
to derive effective models in reduced dimensions such 
as 2D or ID models by excluding polarizations within 
the target layer/chain |8[. With this decomposition and 
within the RPA, the fully screened W can be obtained in 
a two-step procedure as [H 



and 



w = (i- vxir 



w = (i-w x t a r 1 w, 



(i) 



(2) 



where W describes a screened Coulomb interaction ex- 
cluding a specified subset of excitations xb- These exci- 
tations are taken into account when the effective model 
with the interaction W is solved. Alternatively, W 
is obtained from the fully screened W, by rewriting 
Eq. © as 



w =w(i+xbwy 



(3) 



In the present scheme, W corresponds to £7 DMFT and xb 
is a one-center or local target polarization formed at the 
impurity site. 

In practice, the static independent-particle polariza- 
tion formed in the target bands (tb) is calculated using 



Etb r p 

(r,r)=2^^ J ^ U 4(rWr^(r^,(4) 

a/3 qk £ak 

where {V'ak, £ak} are one-body wavefunctions and their 
energies with the wave vector k and the band index a. 
The factor of 2 comes from the spin sum. The band sum- 
mation is performed only over the target bands in the ef- 
fective model. Since the Bloch wavefunctions are related 
to the Wannier functions via the unitary transform as 

1 

7n 



= E ^Si^aW. ( 5 ) 
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the polarization can be recast as 
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X C i R 1 ( r )^R2( r )</ , ofcR3( r ')^R4( 1 ''), (6) 



where m-p, i-l, R1-R4 are the orbital, primitive site, su- 
perlattice site indices respectively and TV indicates the 
total number of superlattice sites. With this expression, 
we specify the target-band polarization formed at the im- 
purity site (the Oth site in R=0) as 

Xo np ( r > r ') =E C 'w ra p^moo( r )0noo(r)0*oo( r ')^oo(r') (7) 



with 



etb 
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nmop 



_ 2 ff*+qrfok ( n Uk) \ * f t(k+q)^ f t(k+qA * rr t(k) 



13 \yoOJ3 J u pO,a 



(8) 



corresponding to the local one-center components of a 
polarization matrix in the Wannier orbital basis. Now, 
by identifying xb i n Eq. © as Xo" 13 an d W as C/ DMFT , we 
write the Dyson equation for the effective interaction as 

W(t,v) = C/ DMFr (r,r') + J> dr"J dr'"C/ DMFT (r,r")x np (r", r") 

xW{r"',r). (9) 

Multiplying this equation by c/4 oo (r)</>„oo(r)0* oo (r') 
x0poo(r') and integrating over r and r', we have 

w, v = u^r + e c c ^- ( io ) 

where we introduce a composite index (/x, v)={(mn), (op)} 
and the matrix element of 0=$V, [/ DMFI } is given by 

O mncp = / dr / dr'^ oo (r)(/ 1 „ O o(r)O(r,r')0* oo (r')</>poo(r / )- 



Thus, Eq. pop is rewritten in a matrix form as 

U DMFT = W(1 + CW)" 1 . (11) 

The equation resembles the unscreening equation ([3]), but 
it is formulated entirely in terms of "local" one-center 
quantities, that can be evaluated straightforwardly, al- 
lowing for a computationally efficient treatment. 

-Application to the Hubbard model. We first apply this 
scheme to the derivation of [/ DMFr for the 2D single-band 
Hubbard model. This is helpful to get insight into the 
behavior of [/ DMFr with respect to changes of the electron 
filling. The Hubbard Hamiltonian reads 

h=-*e c *W-*'E c L c ^-ME niCT+c/ E ni t n 4 ; 

(ij)a ({ij))<r iff i 

where c\ a (cj CT ) creates (annihilates) an electron with spin 
a at site i and n{ a = c\ a a a . t (t') is a transfer integral 
to the (next-)nearest neighbor sites in the (((«, j))) 

sums. U(=8t) and (i represent the onsitc Coulomb re- 
pulsion and chemical potential, respectively. Taking into 
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account the contributions from the charge susceptibility 
only (hence being in accordance with ab initio methods), 
the unscreening equation corresponding to Eq. (fTTjl be- 
comes UM 



U 



DMFT 



u 



w 



I- AW 



(12) 



Here W is a diagonal element of a real-space NxN ma- 
trix W=(l— Uxo) _1 U, TJ a diagonal matrix with ele- 
ments U=U /2 and —A (with A > 0) the diagonal ele- 
ment of the real-space polarization matrix % with ele- 
ments (Xo)ij = Xo(R-i~ R-j)- The latter is obtained by the 
Fourier transform of the reciprocal-space static polariza- 
tion function 



xo(q) 



/k+q - /k 



Y^,/k+q 



(13) 



with £k=— 2i(cos k x +cos k y )~ At' cos k x cos k y — \i and /k 
being the eigenvalue and the Fermi distribution function, 
respectively. 
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FIG. 1: (color online) Filling dependence of [7 DMFr calculated 
(a) with Eq. (fT2")l and (b) with the approximation [Eq. (TH)) ] for 
t' = (red), O.lt (green), 0.2t (blue), and 0.3t (purple). The 
arrows indicate the fillings at which the van Hove singularity 
resides at the chemical potential. 



Figure QJa) shows the filling dependence of £/ DMFT /i 
with various t'. Contrary to a naive expectation, [/ DMFT 
is larger than U . Furthermore, the filling dependence of 
rjDMFT j g no .j. mono tonic and depends on t' . For t'=0, 
rjDMFT a s t ron g peak at half filling where the van 
Hove singularity resides at the Fermi energy. With in- 
creasing t', the peak shifts to higher filling with reduced 
peak height, and another rapid increase emerges at fur- 
ther higher filling. 

These filling and t' dependences of f/ DMFT are 
well understood by the second-order approximation in 
{B n } O: 



JV-l 



U 



DMFT 



u 



U 



XJA" 



(14) 



where B n =[xo)i,i+n is the non-local contribution to the 
polarization. Since the second term of the right hand is 
always positive, the inequality U™ FT >U holds. Figure 
shows the results of [/ DMFT calculated with Eq. (0 
for various fillings and t' . We see in Fig. [TJb) that 
Eq. (fT4"|) well reproduces the overall trend in Fig. [Ha). 

The inequality [7 DMFT >[/ reveals anti-s creemng in- 
duced by non-local polarizations {B n }. This anti- 
screening is intuitively understood as follows: Suppose 
that a test charge electron is put on the impurity site. 
The local polarization screens this electron by creating 
holes at the impurity site. On the other hand, the non- 
local polarizations induce holes or electrons at other sites. 
Then, in the second order process, the induced charges 
create electrons at the impurity site, enhancing the effec- 



tive repulsion. Since 



u 



l + UA 



in Eq. (fT4")l varies smoothly 



with filling [19(, the non-local polarizations {B n } indeed 
dominate the peculiar filling dependence of [/ DMFT . 

In real materials, off-site Coulomb interactions may 
play a role. To see this effect, we have studied [/ DMFT for 
a model with the off-site interaction 1/er with varying e. 
We find that the overall filling dependence of [/ DMFr is 
basically the same as that of the Hubbard model while 
decreasing e (i.e., increasing off-site interaction) causes an 
appreciable reduction of [/ DMFT (not shown). The long- 
range Coulomb interactions connect the onsite polariza- 
tions at different sites and thus bring about the screening 
to the impurity-site interaction. Note that this screening 
works from the zeroth order in {B n }; the approximated 
L/DMFT without the contributions from {B n } indeed be- 
comes smaller than U and has only a weak filling depen- 
dence. 

-Application to SrVO^. We next present ab ini- 
tio results of U DMFT for SrV0 3 . This material is a 
d 1 metal and one of the most benchmarked systems 
within LDA+DMFT (local density approximation plus 
DMFT) 0. On the basis of the DFT band structure, 
we define the target bands by the low-energy t2 g bands 
as was done in Ref. @. We construct three MLWOs per 
V site from the t 2g Bloch states and calculate [7 DMFT for 
these three orbitals. The implementation details and the 
convergence checks are elaborated in Ref. [l8j |. 

Table U compares the values of the onsite intra- and 
inter-orbital Coulomb repulsions (U and U') and Hund's 
rule coupling (J) for the bare (v), cRPA (U cRPA ) [nj, 
U DMFr , and full-RPA (W) interactions. The bare 
Coulomb interactions (~15 eV) are largely screened by 
the high-energy bands, to give J7 cRPA ~3 eV. In the 
present case of SrV03, t/ DMFT turns out to have a value 
similar to [/ cRPA . 

The situation changes drastically, however, when we 
increase the filling n within the rigid-band approxima- 
tion. The left and right panels in Fig. [2] plot U and U', 
respectively, against the filling n. For comparison, we 
also show the results without the non-local polarizations 
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involving the impurity site, i.e., the interaction parame- 
ters calculated without the local one-center and "wing" 
components of the polarization matrix in the Wannicr 
basis ("no- wing" method) [2l|. The result is denoted as 
U no - win s. We see that the filling dependence of U' is 
similar to that of U, except for a constant shift. 

As the filling n increases from 1, [/ DMFT increases 
more rapidly than JJ cRPA . This suggests that the non- 
local anti-screening effect increases more rapidly than the 
screening. Around n=2, {7 DMFT turns to decrease, cross- 
ing [/ cRPA at n^3.5. Finally around the filling end n^5, 
jyDMFT a g a j n increases, as seen in the Hubbard model. 
We see jj no - wins <U rMET at all fillings. This is consis- 
tent with the model analysis: The non-local contribu- 
tions to the screening induce an anti-screening and lead 
to the increase of the onsite interaction. u n °- wln s [ s also 
smaller than JJ cRFA and only weakly depends on the fill- 
ing, consistently with the model analysis where the off- 
site Coulomb interaction induces a screening weakly de- 
pendent on filling. These comparisons clearly show that 
the non-local polarization is the main source of the exotic 
filling dependence of [/ DMFr . 

It becomes now clear that the similar values of {/ DMFr 
and fJ cRPA for SrVC>3 is just a consequence of an approx- 
imate cancellation of the anti-screening by the non-local 
polarizations with the screening by the long-range inter- 
action. In addition, u cRPA ~U DMFT ~U no - wins for SrV0 3 
is partly ascribed to the small filling of the d 1 system 
where the polarization and screening are not large. 

In the previous DMFT studies for the ab initio model, 
rather large values of U compared to [/ cRPA have been 
needed to reproduce the experimental results (e.g., the 



insulating behavior of LaTi03 [24j). Similarly, for the 
2D Hubbard model, the Mott transition takes place at a 
substantially larger U in the single-site DMFT than in 
its cluster extension 
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These aspects are ascribed to 
the intersite correlation effects ignored in the single-site 
DMFT with original [7 cRPA or U . The present scheme 
with {/ DMFT at least partially takes account of the off-site 
effects and will improve the results of the DMFT. The 
vertex corrections ignored in the RPA form have been 
estimated to be small for the conventional cRPA [l[ . For 
the present case, this estimate is left for future studies. 

-Conclusion. We have examined a scheme to evalu- 
ate the effective onsite interaction [7 DMFT for the DMFT. 
Through the analysis based on the Hubbard model, 
we have found unexpectedly an anti-screening effect in- 
duced by non-local polarizations, which competes with 
the screening effects caused by the off-site Coulomb in- 
teraction in real materials. The anti-screening causes a 
non-trivial filling dependence of £/ DMFT and increases the 
effective interaction. Combining the present method with 
DFT, we have indeed shown that [/ DMFT for SrV0 3 ex- 
hibits non-trivial filling dependence if the chemical po- 
tential is varied. 

-Acknowledgments. We would like to thank Takashi 



TABLE I: Onsite bare (v), cRPA (U cRPA ), present-scheme 
(U DMFr ), and full- RPA (W) interaction parameters calcu- 
lated for SrVC>3. The unit of energy is eV. The method was 
implemented in two codes, Tokyo Ab initio Program Pack- 
age [x| (lef t values) and the Vienna Ab initio Simulation 
Package [23| ] (right ones), which yield almost identical values 
for U cRP . Otherwise, the latter values are generally 5-10 % 
larger than those of the former, since the exact shape of the 
orbitals is used in VASP. 
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FIG. 2: (color online) Filling dependence of intra-orbital 
(left) and inter-orbital (right) screened Coulomb repulsion 
of SrVOa evaluated within full RPA, cRPA, present scheme 
(U DMFT ), and "no-wing" methods, which are calculated with 
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Supplemental Materials 



S.l DERIVATION OF EQS. (12) AND (14) 

An RPA fully-screened interaction W may be ex- 
pressed as 

W = e- 1 U, e=l-Ux . (S.l) 

Here X = I W, U, % , e] are 2Nx2N matrices decom- 
posed into their spin channels according to 



With this decomposition, U and Xo are written as 



U = 



u 

U 



and Xo = 



hx o 

3X0 



(S.3) 



respectively, where U is a diagonal matrix with elements 
U and Xo is a real-space polarization matrix. In the Hub- 
bard model, only the onsite components W = [W, 
are relevant, which are given by 



UJ ii 



w = 



e" x U 

TT 



t4- 4-4- 



(S.4) 



According to Eqs. (5.9) and (5.10) in Ref. [Sl[, the in- 
verse dielectric matrix in the ft spin channel is (1 — 

(J-UXo) 2 ) 1 so that we obtain 



W 



'\ V Xof 



U 



This equation is also written as 



W = 



u + iu( Xc -x s )u 



(S.5) 



(S.6) 



with Xc (x.s) being the charge (spin) susceptibility given 

by X c =(l-ix U)- 1 i Xo [Xs=(l+5XoU)- 1 |xo]- In 
line with ab initio methods, which only take charge fluc- 
tuations into account, we consider the term related to 
Xc om yj the resulting expression for W is 



W-- 



U + iu Xc U 



u- 



u 



1 - ~ 2 XoV 



U + U(l- Xo U) 'xoU 



U+(l + U Xo + (U Xo ) 
U+(l-U Xo ) _1 U 

u + w, 



)U 



(S.7) 



where U=U/2 and W=[(l-tJx ) U]«- 

We now decompose the total polarization Xo m to the 
two parts, 



Xo 



-A0\ , f B T 

o) and Xo= (b xo' 



(S.8) 



where B=(i?i,i32, ••• ,-Btv-i) T and x'o 1S an (TV — 
l)x(TV — 1) matrix. Then, replacing Xo with x'o m 
Eqs. (fSH -7). we obtain 



jjFMFT _ jj , jjTMFT 



(S.2) with 



U 



DMFT 



(1-UXo) U 



(S.9) 



(S.10) 



ii 



The above derivation of {7 DMFT is based on the screening 
approach of Eq. (2). On the other hand, [7 DMFT can also 
be obtained in the unscreening approach of Eq. (3) as 



V 



DMFT 



W(1 + X *W) 



w 



I- AW 



(S.ll) 



Eqs. (fST9|> and (jS.lip give Eq. (12) in the main text. 
Again using Eqs. (5.9) and (5.10) in Ref. |l|, Eq. (ISlOl) 
is further recast into 



U 



DMFT 



1 - U 2 B T {1 - U X 



= U. (S.12) 



Hence, up to the second order in {B n }, we obtain 



jjYMFT ^ jj 



1 + 



1 + UA 



(S.13) 



which is equivalent to Eq. (14). 



S.2 IMPLEMENTATION DETAILS OF U DMFr IN 
THE PLANE- WAVE BASIS-SET CODE AND 
COMPUTATIONAL RESULTS 



Here, we describe implementation details for the ab 
initio U DMFr calculations. The calculation is per- 
formed with the norm-conserving pscudopotcntial and 
plane-wave basis set and the projector augmented wave 
method, respectively [S2j, S3]. In the plane- wave basis-set 
calculation, two different cutoffs for the plane waves are 
conventionally used; the low-momentum cutoff g'°^ for 
the polarization function and the high-momentum cutoff 
gcif h f° r orbitals. In general, the structure of the polar- 
ization function in real space is smooth compared to that 
of the wavefunction, so we can employ the smaller cut- 
off and it considerably reduces the computational cost. 
In the U DMFr calculation in Eq. (11) in the main text, 
however, we should be careful about the use of the two 
different cutoffs. 
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The Dyson equation Eq. (10) is written in the momen- 
tum space with the double Fourier transform |S4j as 



^ glg2 ^^E^^T g4 ^ g2 (l&l<gc°uT)(S.14) 



g3g4 

W _r?DMFr r rsr low <|cr |<B- Mgh ' 

vy glg2~ t7 glg2 ~ u gl°glg2 \&cut — l&il— Bcut 



), (S.15) 



where gi-g4 are reciprocal wave vectors associated with 
the superlattice jsi[ and v s ~4:ir/\g\ 2 is the Fourier trans- 
form of the bare Coulomb interaction v. In Eq. (|S.15jl we 
have used the fact that x™f vanishes outside gj.°™. 

Recognizing this aspect, we define the low- and high- 
momentum contributions for W^ Ul defined in Eq. (10), 
as 

-. low 

^-T7E(^^ol eiSr ^^o) W ss'(^ol e " iS ' r '^oXS.16) 



gg' 

high 



=7/X>™°°l e * Sr V>"00 ) w s ^°oo| e - l&r V P oo) (S . 17) 



Here, y is the c rystal volume and W^Wj^+W^. 
The sum in Eq. (|S.14[) is taken for the reciprocal vector 
within g °^, while the sum in Eq. (|S. 15[) runs over the 



reciprocal vector for gc°7<|g|<g"f. Similarly, [7°^ is 
written as the sum of U^f 1 ^ and Uf^-^ . Insert- 
ing Eq. (IS.14j) into Eq. (|S.16j) with the double Fourier 
transform of X mp , we obtain 



Ti^low rrDMFT-low 

y * liv ^ [III 



or in the matrix form 



•^ylow -jj-lJMh'l -low _|_ -^jJJlVlh'l-low^-^y-low 

Since u DMFT =U mlFT - low +U Dwr - hi g h , after some manip- 
ulations, we obtain 



tDMFT-Iow 



tDMFT-Iow/ 



U 



DMFT 



W'™(1 + CW 



low \ 



yhigh 



(S.20) 



with V high (=w hi s h =U DIV]Fr - hi s h ) being the matrix of v 
at high momenta Eq. f|S . 1 7[) . In the actual calculation, 
this expression is used. 

As a note on the numerical calculation, we remark 
some details for calculating the polarization function in a 
metallic system. The target-band polarization Xg b (r, r') 
in Eq. (4) in the main text is given in the momentum 
space with the double Fourier transform as 



' (q) - 2 E E (^k| e - t(q+G) -i^ 



AGG 



k a/3 



: (^/3k 



+ql 



k+q/ 

' (q+G,) - r '|Vak>. (S.21) 



Here, G is a reciprocal lattice vector for the primitive 
lattice and q is a wave vector in the first Brillouin zone. 



{"0ak}, {eak}, and {/„k} are the Bloch states, their ener- 
gies, and occupancies, respectively, and the band summa- 
tion runs over the target bands only. In the calculation 
of Xgc (q) °f the- metallic system, the k integral on the 
right hand must be performed carefully, because the ex- 
pression includes a numerical instability due to the Lind- 
hard part. To avoid the instability, we use the Wannier 
interpolation scheme }S6l |: we interpolate the original k- 
point data (of about lOx lOx 10) for the eigenvalues {e Q k} 
and interstate matrix elements {(i/'/3k+q|e^ q+G ^ r |'!/ , ak)}, 
to obtain the data on a denser k grid (about 30x30x30). 
After such an interpolation, the k integration is per- 
formed with the generalized tetrahedron method [S7[ to 
obtain both, real and imaginary parts of XgG'W- 

We also need a careful treatment of poles at e^jk+q = 
e Q k in Eq. ()S.21[) , for which we rewrite 



fpk+q — foik 



Cak 



e /3kH 



e Q k 



(S.22) 



Based on the central-difference approximation of the 
Fermi-Dirac function with The Fermi level ep- Switch- 
ing to the 6 function in Eq. (|S.22p is performed in the 
threshold |e^k+q — e a k|<0.06 eV and the 6 function is 
treated with a smearing factor of 0.03 eV. With the re- 
sulting targ et-band polarization \ tb and the rest polariza- 
tion x r [S8I ] , the fully screened RPA Coulomb interaction 
W low in Eq. (jS^O]) is calculated, where the W GG >(q)m- 
teraction at q — > limit is treated following Ref. [Slj . 
The same treatment is applied to the evaluation of the 
Wannier matrix elements of (7^ in Eq. (8) [S9| | . With 
all these treatments, the present XJ DMFT calculation en- 
sures the accuracy within several percent. 

If not otherwise noted, the density-functional the- 
ory calculations for SrV03 were perfo rmed with 
Tokyo Ab initio Program Package [S10(, which is 
based on the pscudopotcntial plus plane-wave frame- 
work. The exchange-correlation functional is calcu- 
lated within the generalized-gradient approximatio n with 



Perdew-Burke-Ernzcrhof (PBE) parameterization jSll 



and t he T roullier-Martins norm-conserving pseudopoten- 
tials S12 1 in the Kleinman-Bylander representation |S13| 
is adopted. In the present calculation for Fig. 2 in the 
main text, the cutoff energies for wavefunctions and po- 
larization functions are set to 49 Ry and 25 Ry, respec- 
tively, and we employ llxllxllk points. The Brillouin- 
zone integrals are evaluated using the generalized tetra- 
hedron method [s3] after interpolation to a 33 x 33 x 33 k 
mesh. 

Where noted, additional calculations were performed 
using the Vienna Ab initio Simulation Package (VASP), 
using projector augmented waves and the local density 
approximation. The plane wave cutoff energies for the or- 
bitals and response functions were set to 414 eV (30 Ry) 
and 250 eV (18 Ry), respectively. Extrapolation to a high 
energy cutoff (500 eV) was performed using Eq. (|S.20[) . 
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In VASP no intermediate extrapolation to a denser k- 
point grid was performed. Instead, in Eq. (|S.21[) . the 
Fermi occupancy function /(e) was replaced b y a M eth- 
fessel Paxton smearing function with cr = 0.1 [S14j . and 
consistent with metallic screening Woo(q — > 0) was set 
to 0. 

Figure S Q] shows our calculated band structure of 
SrVC>3 (a) and the density of states for the t2 g bands 
(b). The arrows in the panel (b) indicate the Fermi lev- 
els for the fillings n=1.0 to 5.0 with the interval 0.5. We 
see that the van Hove singularity nearly corresponds to 
the Fermi level at the filling n = 4.0. 




Energy (eV) 

Fig. S 1: (color online) (a) Calculated electronic band struc- 
ture of SrVC>3. The interpolated band dispersions for the t2 g 
bands are depicted as blue dashed lines, which cross the Fermi 
level, (b) Calculated density of states for the t2 9 bands. Black 
arrows indicate the Fermi level for the filling n = 1.0-5.0 from 
left to right for the values shown in Table S IIVI 



We show in Table S U and S |TT] the convergence behav- 
ior of U DMFr calculated for SrVC>3 against the sampling k 
points using the Tokyo Ab initio Program Package. The 
table lists the values for the onsite intra- and inter-orbital 
Coulomb repulsions (U and V) and Hund's rule coupling 
(J). The usual constrained random-phase-approximation 
(cRPA) (U cRPA ) [sl and full-RPA (W) results are also 
shown for comparison. We see that the results are al- 
most converged at 6x6x6 or 7x7x7 k-point samplings. 
Despite a less sophisticated interpolation procedure the 
results using the Vienna Ab initio Simulation Package 
(VASP) show a very similar convergence behavior. Again 
the error is reduced to few percent at 7x7x7 k-points, 
although a sizeable scattering prevails in both codes. 



Table S IIIII shows the convergence behavior against 
the cutoff momentum g|.°™ for the polarization function. 
We see that the convergence is attained around gJ.°Y~25 
Ry. Finally, Table S |TV] lists the interaction parameters 
calculated at the fillings n—1. 0-5.0, which are used for 
the plot in Fig. 2 in the main text. In this table, we 
add the "no-wing" data (U n °- win s). For the definition of 
U no " win s, see the main text. 
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Table S I: Convergence behavior of U cRPA , XJ DMFT , and W to the sampling k points of SrVC>3 for the Tokyo Ab initio Program 
Package. The cutoff energy for polarization function is 25 Ry. 



17 VP J u VP J u VP T 



5x5x5 3.40 2.34 0.47 3.48 2.41 0.41 0.93 0.23 0.33 

6x6x6 3.50 2.45 0.47 3.44 2.37 0.47 0.98 0.25 0.33 

7x7x7 3.42 2.37 0.47 3.37 2.30 0.47 0.97 0.25 0.33 

8x8x8 3.32 2.27 0.47 3.26 2.20 0.48 0.96 0.24 0.33 

9x9x9 3.27 2.22 0.47 3.22 2.16 0.48 0.97 0.25 0.33 

10x10x10 3.44 2.38 0.47 3.39 2.33 0.47 0.98 0.25 0.33 

llxllxll 3.39 2.34 0.47 3.33 2.27 0.47 0.97 0.25 0.33 



Table S II: Convergence behavior of U cRPA , U DMFT , and W to the sampling k points of SrVC>3 for the Vienna Ab initio 
Simulation Package. 



U cHFA U UMbT W 

U VP J U U' J U VP J 

3x3x3 3.45 2.43 0.50 6.38 5.38 0.48 1.02 0.23 0.38 

4x4x4 3.31 2.30 0.49 5.25 4.26 0.47 1.00 0.22 0.38 

5x5x5 3.31 2.30 0.49 3.94 2.95 0.47 1.07 0.26 0.39 

6x6x6 3.35 2.34 0.49 3.50 2.51 0.47 1.11 0.29 0.39 

7x7x7 3.38 2.36 0.49 3.51 2.53 0.47 1.17 0.34 0.40 

8x8x8 3.36 2.35 0.49 3.46 2.47 0.47 1.12 0.30 0.39 

9x9x9 - - - 3.42 2.43 0.47 1.10 0.29 0.39 

10x10x10 - - - 3.42 2.43 0.47 1.11 0.30 0.39 

llxllxll - - - 3.48 2.49 0.47 1.14 0.31 0.39 



Table S III: Convergence behavior of U cRPA , U DMFT , and W to the cutoff energy for polarization function gj.°™ f° r the Tokyo 
Ab initio Program Package. The sampling k points are fixed at 7x7x7 and, in the interpolation of the polarization calculation, 
the 21x21x21 fc-grid is employed. 



jjcRPA jjDMFT w 

u VP J u VP J u u' j 

10 Ry 3.48 2.37 0.51 3.38 2.28 0.51 1.22 0.26 0.45 

15 Ry 3.48 2.39 0.49 3.39 2.30 0.49 1.13 0.27 0.39 

20 Ry 3.44 2.38 0.48 3.37 2.30 0.48 1.04 0.26 0.36 

25 Ry 3.42 2.37 0.47 3.37 2.30 0.47 0.97 0.25 0.33 

30 Ry 3.41 2.36 0.47 3.37 2.30 0.47 0.94 0.24 0.32 

35 Ry 3.40 2.36 0.47 3.37 2.30 0.47 0.91 0.24 0.31 
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Table S IV: Our calculated U cRPA , U DMFr , and W at fillings n=1.0-5.0 (Tokyo Ab initio Program Package). These data are 
used in Fig. 2 in the main text. The u no " wmg data are also listed. For the definition of u no_wms , see the main text. 



■jjcRPA 


■yDMFT 


jj-no-wing 


w 


U U' J 


U U' J 


U U' J 


U U' J 



n=1.0 3.39 2.34 0.47 3.33 2.27 0.47 3.30 2.24 0.47 0.97 0.25 0.33 

n=1.5 3.47 2.41 0.47 4.01 2.93 0.48 3.36 2.29 0.48 0.80 0.16 0.29 

n=2.0 3.65 2.59 0.46 4.74 3.63 0.47 3.41 2.34 0.47 0.68 0.11 0.26 

n=2.5 3.72 2.65 0.46 4.58 3.48 0.47 3.23 2.16 0.47 0.59 0.07 0.24 

n=3.0 3.83 2.75 0.45 4.33 3.23 0.46 3.14 2.07 0.46 0.53 0.06 0.22 

n=3.5 3.89 2.81 0.45 3.85 2.76 0.45 3.01 1.96 0.45 0.49 0.04 0.20 

n=4.0 3.93 2.85 0.44 3.39 2.32 0.44 3.02 1.96 0.44 0.47 0.04 0.20 

n=4.5 3.98 2.90 0.44 3.05 2.00 0.43 2.94 1.90 0.43 0.50 0.05 0.20 

n=5.0 4.06 2.97 0.43 3.58 2.50 0.43 2.75 1.71 0.42 0.62 0.08 0.24 



